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' p • Abstract 

We report the theoretical evaluations of the static scalar polarizability of the 133 Cs ground state 
and of the black body radiation shift induced on the transition frequency between the two hyper- 
fine levels with tuf = 0. This shift is of fundamental importance in the evaluation of the accuracy 

•l-H . 

of the primary frequency standards based on atomic fountains and employed in the realization 

i-£h ' of the SI second in the International Atomic Time (TAI) scale at the level of 1 x 10 -15 . 

Our computed value for the polarizability is Qo = (6.600 ±0.016) x 10~ 39 Cm 2 /V in agreement 

at the level of 1 x 10~ 3 with recent theoretical and experimental values. As regards the black body 

radiation shift we find for the relative hyperfine transition frequency = (—1.49 ± 0.07) x 10~ 14 

at T = 300 K in agreement with frequency measurements reported by our group and by Bauch 

and Schroder [Phys. Rev. Lett. 78, 622, (1997)]. This value is lower by 2 x 10~ 15 than that 
CO 

obtained with measurements based on the dc Stark shift and with the value commonly accepted 



up to now. 



PACS number(s): 32.60. +i, 32.10.Dk, 06.30.Ft 
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I. Introduction 



The evaluation of the accuracy of the cesium primary frequency standards requires to take into 
account, among others, the shift induced on the clock transition \6 2 Si/2,F = 4,niF = 0) — > 
\6 2 Si/2,F = 3,mp = 0) by the black body radiation (BBR). The definition of the second in the 
International System of Units (SI) is in fact based on the above transition observed in absence 
of any perturbation and then, in particular, at T = K [Q, where T is the temperature of the 
environment. 

Gallagher and Cook [2] firstly considered the effects of BBR on the Rydberg states of the 
atoms while Itano et al. j3j focused their attention onto the atomic frequency standards. They 
reported the following expression for the relative frequency shift of the 133 Cs clock transition 
(i/ = 9192631770 Hz): 

with = -1.69(4) x 1CT 14 and e = 1.4 x 1(T 2 when T = 300 K. This expression has been 
adopted till now to correct all the primary frequency standards for the BBR shift at an accuracy 
level of 10 -15 or even better [I]. 

The BBR induces both non resonant ac Stark and Zeeman shifts of the atomic transitions. 
For the ground state 2 Si/2 of the alkali-metal atoms the Zeeman shift is of the order of 10~ [2] 
and is not relevant for the accuracy estimation of the present frequency standards. 

The direct experimental method to evaluate the BBR shift is based on (fTJ): the frequency shift 
Av is measured versus the black body temperature T and the experimental points are fitted with 
(^Q) whose temperature dependence follows directly from the Planck radiation law. This method 
has been used by Bauch and Schroeder [Sj employing a thermal cesium beam and in our group 
using an atomic fountain 0. 

In order to increase the experimental resolution, the BBR shift may be also evaluated indi- 
rectly by measuring the dc Stark shift coefficient k defined as: 

Av = kE 2 (2) 

where Av is the frequency shift induced by the static electric field E. Measurements of k have 
been performed in this way by Haun and Zacharias [7] and by Mowat using a thermal beam 
and the Ramsey interrogation technique, and by Simon et al. jHJ using an atomic fountain. The 
BBR shift is then obtained from © assuming E: 

E 2 = (£ 2 (t)) = ^l(^) 4 (3) 
e c \T J 

where (£ 2 (t)} is the mean-squared electric field of the black body radiation, a = 5.670400(40) x 
10 -8 W/m 2 K 4 is the Stefan-Boltzmann constant jlOj . eq is the vacuum permittivity and c is 
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the speed of light in vacuum. In (j3J) the first identity is valid at a level of a few percent and 
can be corrected as reported in [3], while the second identity is obtained directly by the Planck 
radiation law (see next section). 

The measurement of the scalar static polarizability oq is related to the above method through 
the expression: 

W = -\a G E 2 (4) 

where W is the potential energy of a neutral atom in a static electric field E of moderate strength. 
Due to the very low intensity of the black body field, the hyperpolarizability effect is negligible 
[II] and will not be considered throughout the paper. 

The polarizability is obtained from the force acting on the atom when submitted to a gradient 
of the electric field: 

F = -VW = a EVE (5) 

This force, in turn, can be measured through the deflection of thermal beams |12| I13|. through 
E-H gradient balances |14|I15| or through the time of flight of cooled atoms in a fountain apparatus 
|16| I17j. In principle, the BBR shift can be also obtained from Q and (jU as follows: 

A„ = 1 -(W P - W a ) = -JL [ aom - a (\a})) (£ 2 (t)} (6) 

where \/3) = \6 2 Si/2,F = 4, mp = 0) and |a) = \Q 2 Si/2,F = 3,mp = 0) in the Cs case; h is the 
Planck constant. 

The indirect experimental technique allows a higher resolution, being possible to submit the 
atoms to dc fields of several MV/m, if compared to the black body field W (£ 2 (t)) = 832.2 V/m 
at T = 300 K used in the direct method. Nevertheless, it is based on the conceptual assumption 
E 2 = (£ 2 (t)) and may not be free from spurious effects as we shall discuss further on. 

^From the theoretical point of view, the evaluation of the BBR shift requires the computation 
of the scalar polarizabilities of the two states defining the clock transition or, in other words, of 
the non resonant ac Stark shift. This is typically performed through perturbation techniques 
|18l I19j . More precisely, according to ©, we have to calcolate the differential polarizability 
between the two ground state hyperfine levels. As pointed out by Feichtner et al. j5U], this 
requires to consider in the theoretical analysis a basis of modified eigenfunctions which account 
for the hyperfine interaction of the ground and P states with other states of the same quantum 
numbers F and mp. A review of the theoretical works leading to the evaluation of the ground 
state 133 Cs polarizability can be found in JH] and in |17j . 

Due to the basic importance of the BBR shift correction to evaluate the present primary 
frequency standard accuracy at the 10 -15 level and, in the next future at the 10 -16 level, it is 
highly desirable to re-examine the relation considering that: 
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(i) new and more precise experimental data are now available for the electric dipole moments 
and for the transition frequencies of 133 Cs involved in the theoretical computation of ao; 

(ii) recent theoretical evaluations of ao l ea, d to a significantly lower value than that com- 
monly assumed up to now [HI IE] ; 

(iii) the two most precise measurements El of the electrical polar izability do not agree to 
each other. 

In this paper we report an ab-initio computation of the BBR shift for the clock transition which 
turns out lower by 2 x 10~ 15 at 300 K than the commonly adopted value of (^Q). A critical 
comparison with the theoretical and experimental values reported in the literature will be also 
given, showing the present difficulty to asses the accuracy at the 10 -16 level for atomic frequency 
standards operating at T ~ 300 K. 



II. Basic Theory 

The non-resonant ac Stark shift AWa of the energy of state \a) induced by the level \k) is given 
by [S]: 

AW (k) = \dakl_ [ + °° v a -v k g 2r v)dv (7) 

where d a k = e(a|e • r\k) is the electric dipole moment of the transition \k) — ► \a), e the electron 
charge, e the polarization vector of the electric field, r the atomic position vector, v a = W a /h 
and Vk = Wk/h the energy levels of the states \a>) and \k) in frequency units and £„(y) the mean- 
squared value of the electric field at frequency v. When \a) represents a ground state, we have 
v a — v k = ~ v ak = ~ c /^ak being v a ^ and the frequency and the wavelength of the atomic 
transition \k) — > \a). Figure 1 summarizes the above definitions. The electric field per frequency 
unit £ u {v) is related to the spectral energy density E u (i , ,T) by the well known relation: 

£ 2 v (v)dv = —E u (u, T)dv (8) 

In our case the energy density E v {u,T)dv is given by the Planck radiation law: 

8i:hv 3 dv 

E v {y, T)dv = c 3 ehv/kBT _ 1 (9) 

where ks is the Boltzmann constant. 

The maximum of E u (v,T) is reached when (Wien law): 

v = iw = 2.821^T (10) 

At T = 300 K, relation gives v max — 17.6 THz. In the Cs case, the minimum value of 
v a k corresponds to the Di optical transition (335 THz) so that v max <C v a k for all \k) levels 
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Figure 1: Atomic levels considered in (J7J). 



coupled to the state |a) and it is possible to omit the v 2 term in the denominator of ® avoiding 
the singularity of the integrand at v = u a k. We shall consider further on the effect of this 
approximation. The basic relation Q can then be written as: 

AW (k) = * \^L^ { t) ) (ii) 
h v a k 

where the mean-squared electric field (£ 2 (t)) is given by 

(£ 2 (t)) / £ 2 v {v)dv = — (12) 

The relation (|12jl. obtained through (jHJ) and @, is also known as the Stefan-Boltzmann law; at 
T=300 K it gives (£ 2 (t)} = (832. 2) 2 (V/m) 2 , as already reported in the previous section. 

Summing over all the excited states \k) coupled to the level |a), we obtain the total Stark 
shift for the state \a): 

A^« = -££^V(t)> (13) 
Moreover, taking Q and Q) into account, it is possible to write (|T3|) also in the following form: 

AW a = ~a (£ 2 (t)) (14) 

where 

2 \d a k\ K \ 
«o = 7> (15) 

is the scalar polarizability of the state |a) which may be measured experimentally. 

More precisely, the expression (|15|) takes into account the contribution to the polarizability 
from the valence excited states. In principle, the total polarizability includes also the contribution 
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from the core excited states |22| : however, this is not important for the clock transition shift since 
it is common to both the ground state sublevels. 

The BBR shift Av of the clock transition \(3) — > \a) we are looking for is: 

Au = ^(AWg - AW a ) (16) 

Introducing relations (|12[) and (|13|) in (|16|) we finally obtain: 



h 2 e c\T J \^ vpk v a k 

wi/3))-«o(H)} 



(17) 



he 

We conclude this theoretical review observing that the dc Stark coefficient k defined in (J2J) can 
be expressed through the polarizabilities of the ground state levels (see also (JBJ) as: 

k = ~{aom) ~ «o(|a»} = (18) 

where a\o is the scalar differential polarizability introduced by Sandars [22]. In the experimental 
measurement of k performed with a dc electric field a tensorial component is also present, which 
has to be subtracted in order to take the isotropy of the black-body radiation into account [HIE]. 

III. Numerical Evaluation 

To evaluate numerically the BBR shift given by (|17j) we need the matrix elements of the operator 
e • r; writing the atomic position vector r as an irreducible rank 1 tensor operator we have: 

e • r = r cos -=r\ sin 9e~ tlp + — =r_i sin 6e llf (19) 

where r is expressed in terms of its spherical components 7*0, r_ 1 and t\\ (co-latitude) and ip 
(azimuth) are the polar angles defining the electric field direction with respect to the quantization 
axis (z axis). The well-known Wigner-Eckart [21] theorem allows us to simplify the generic matrix 
element: 

( 6Si/ 2 ;F,m F \er q \nPjr,F',m' F ) = Q(F,m F ; F' ,m' F ; J, J'; q){ 6S 1/2 \ |er| \n'Pj>) (20) 

where the double bars indicate the dipole reduced matrix element and q is the index labelling 
the component of r. In ()2Uj) we defined: 

Q(F,m F ;F',m F ,;J,j' ]q ) = ( _ 1) F-m F +F'+.7+/+y (2F + 1){2F , + 1} 

F 1 F' \ $ F 1 F' ) C2\.) 
-ra F q m' F I I J' I J 



where the coefficients in round and in curly brackets are the 3-j and 6-j symbols respectively. 
The primed quantum numbers refer to the excited states \n'Pj*) coupled to the ground state by 
allowed electric dipole transitions. 

The polarizability of each of the two ground state levels can then be written as: 

" f da Y. l(6S " 2lM "' Pj ' )|2 



a (\F,m F )) 



Anh 



n/>6 



2 ^n>Pj,;F>) 



Q(F,m F ; F', m' F ; J, J'; q = 0) cos 9 



1 



F 

7' ™/ 



y=Q(F,m F ; F', m' F ; J, J'; q = 1) sin( 

v2 



-up 



+ - r Q(F,m F ;F',m , F ;J,J';q = -l)sm9e i{p 
v2 



(22) 



In (|22|) we perform an integration over the solid angle to take the isotropic nature of the BBR 
into account, consequently the integral is normalized to 4-7r; v(6Si/2 — ► n'Pj>; F') is the frequency 
of the transition indicated in parentheses. According to selection rules, for fixed values of the 
quantum numbers only one of the three Q terms is not zero. 

Actually, as Feichtner et al. [20] pointed out, the S\/2 state and the P states wave functions 
are perturbed by the effect of their hyperfine interaction with other states with the same F 
and m F but with different principal quantum number n. Therefore, the dipole matrix elements 
should not be calculated between the unperturbed wave functions but between the perturbed 
ones defined as: 



(65 1/2 ; F\ = (6S 1/2 ; F\ + ^2 a n i F {nS 1/2 ; F\ 

n=7 

oo 

(6Pjr,F , \ = (6Pjr,F'\+Y,Kj>F>(nPj>;F'\ 



n=l 

oo 



(23) 



(7Pjr,F'\ = {lPjr,F'\+Y,CnJ'F'{nPjr,F'\ 



n=6 
n^7 



The coefficients a n i F , b n j F and c n j F , given in [2U], are here reported for convenience. 



dnJF 
b n JF 



V>ns(0) 
V>6S(0) 

Ues(0)y 



33/2-F(F + l) 
8 



v(nS 1/2 -> 651/2) 



63/4 + J(J + 1)-F(F + 1) 
32J(J + 1) 



'A) 



v{nPj -► 6Pj) 



(24) 



C n JF 



(V-6S(0)) 2 



63/4 + J( J + 1) - F(F + 1) 
32J(J + 1) 



'A) 



v{nPj -> 7Pj) 



and their numerical values are reported in Tables 1, 2 and 3. 
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Table 1: Coefficients a n jF (in units of 10 6 ) for perturbed cesium 6S wave function. 



n J\F 


4 


3 


7 1/2 


-3.37 


+4.33 


8 1/2 


-1.63 


+2.08 


9 1/2 


-1.04 


+1.34 



Table 2: Coefficients b n jF (in units of 10 6 ) for perturbed cesium 6P wave function. 



n 


J\F 


5 


4 


3 


2 


7 


1/2 


\ 


-1.966 


+2.529 


\ 


8 


1/2 


\ 


-0.904 


+1.16 


\ 


9 


1/2 


\ 


-0.567 


+0.729 


\ 


7 


3/2 


-1.223 


-0.0583 


+0.874 


+1.573 


8 


3/2 


-0.560 


-0.026 


+0.400 


+0.721 


9 


3/2 


-0.351 


-0.017 


+0.250 


+0.452 



Table 3: Coefficients c n jF (in units of 10 6 ) for perturbed cesium 7P wave function 



n 


J\F 


5 


4 


3 


2 


6 


1/2 


\ 


+1.966 


-2.529 


\ 


8 


1/2 


\ 


-1.539 


+1.980 


\ 


9 


1/2 


\ 


-0.741 


+0.953 


\ 


6 


3/2 


+1.223 


+0.0583 


-0.874 


-1.573 


8 


3/2 


-0.957 


-0.0455 


+0.684 


+1.230 


9 


3/2 


-0.455 


-0.0217 


+0.325 


+0.585 



These coefficients are very small (of the order of 10~ 6 ) and their contribution to the absolute 
value of the scalar polar izability is completely negligible. However, they must be taken into 
account in order to obtain the correct value of the differential polarizability we are interested in. 

The ground state wave function is more sensitive to the hyperfine interaction so in this 
calculation the corresponding sum in (J23|) has been extended from n = 7 to n = 9; the coefficients 
b n jF and c u jf essentially involve higher order corrections and we considered only the first term 



8 



in the respective perturbative sums. 

As an example, we report in the following equation the expression of the polar izability of the 
level \6Si/ 2 ;F = A,mp = 0) due only to the interaction with the excited level I6P1/2): 



Am , , , 

t 1 ,m F 



Q(4, 0; F', m' F ; 1/2, 1/2; q = 0) cos 9 



1 -1 ■ 2 

^=Q(4, 0; F , m' F ; 1/2, 1/2; q = 1) sin 6e~ lip + ^Q(4, 0; F', m^; 1/2, 1/2; g = -1) sin 0e^ 



1 



( ^S 1/2 \\er\\6P 1/2 ) + ^a n i A { nS 1/2 \ \er\ |6P 1/2 ) + 



n=7 



K65 1/2 ^6P 1/2 ;F'; 
6 7 i F ,( eSi/alHlTPi/a) +^a n i 4 ^ 7 i F /( nSx/alH^P^) 

n=7 

(25) 

For both the ground-state hyperfine levels, we calculated similar contributions due to the excited 
states |6P 3 / 2 ), |7Pi/ 2 ), |7P 3 / 2 ), |8Pi /2 ) and |8P 3/2 ). In the case of |6P 1/2 ) and |6P 3/2 ) (Di and 
D2 line, respectively), we also included their hyperfine structure |25j . 

To evaluate numerically the absolute and the differential polarizability we need to know the 
reduced dipole matrix elements and the frequencies of the transitions involved in this calculation. 
Several experimental ^3 EHJ [23 EH] and theoretical [Ml EH EU E2] works have been devoted to 
the determination of the dipole matrix elements. In particular, the measurement of the excited 
states I6P1/2) an d I6P3/2) lifetimes provides a direct method to know the reduced matrix elements 
through the well-known expression of the spontaneous emission: 

being tji the lifetime of the transition \k) — > \a) . For Di and D2 lines, performing an average of 
the lifetimes reported in Table 2 of [Hj we obtain: r 1/2 = (34.86±0.05) ns and r 3/2 = (30.44±0.04) 
ns. However, due to the dispersion of these values we considered an error bar that is twice the 
calculated one. For the reduced dipole moments (6S1/2 1 \er\ \ nPj) (n=7,8), (7S±/ 2 \\er\\7Pj) and 
(6Pj\\er\\7Si/ 2 ) we perform an average over the theoretical and experimental values reported in 
Table 6 of [3U]. For the matrix elements (7<S , i/ 2 ||er||8Pj) we use the values of [21] to whom we 
attribute an uncertainty of the 1%. The signs of the matrix elements are chosen according to 
the Feichtner et al. paper. All the values of the dipole matrix elements used in this numerical 
evaluation are summarized in Table 4. 

In Table 5, we report the values of the atomic transitions frequencies we use. For the frequen- 
cies of Di and D 2 lines we use the most accurate values presently available [2S]. For the other 
transitions we use the values reported in the Basic Atomic Spectroscopic Database of NIST. 
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Table 4: Reduced dipole matrix elements of the atomic transitions used in this calculation. 



Dipole Matrix Element 


Value in 10 29 C • m 


(GS 1/2 \\er\\6P 1/2 } 


-3.8174(56) a 


(6<Si/ 2 |M|6P 3/2 ) 


-5.3729(70) a 


(65 1/2 ||er||7Pi/2) 


0.237(4) b 


(6Si/ 2 ||er||7P 3/2 ) 


0.491(3) b 


(6Pi /2 ||er||75i /2 ) 


3.59(l) b 


(6P 3 / 2 ||er||75i/ 2 ) 


5.487(8) b 


(75 1/2 ||er||7P 1/2 ) 


-8.71(ll) b 


(7S , 1 / 2 ||er||7P 3 / 2 ) 


-12.11(l) b 


(6P 1/2 ||er||8S 1/2 ) 


-0.85(9) d 


(6P 3 / 2 ||er||85i/2) 


-1.23(12) d 


(65 1/2 ||er||8P 1/2 ) 


-0.066(2) b 


(651/2 1| er||8P 3/2 ) 


-0.183(2) b 


(75 1/2 ||er||8P 1/2 ) 


0.776(l) c 


<7S 1/2 ||er||8i> 3/2 > 


1.374(14) c 


(7P 1/2 ||er||8S 1/2 ) 


7.41(70) d 


(7P 3/2 ||er||85 1/2 ) 


11.5(l) d 


(6P 1/2 ||er||9S 1/2 ) 


0.46(5) d 


(6P 3/2 ||er||9S'i /2 ) 


0.65(6) d 


(7P 1/2 ||er||9S 1/2 ) 


-1.67(17) d 


<7P 3/2 ||er||9S 1/2 ) 


-2.2(2) d 


a see text 



b from Ref. [H| 
c from Ref. [32] 
d from Ref. [H| 
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Table 5: Frequencies of the atomic transitions used in this calculation. 



Transition 



Frequency (THz) 



651/2 


F 




4>-> 


6^1/2;^' 


= 3) 


335 


1113702 


65*1/2 


F 


= 


4)- 


6Pi /2 ;F' 


= 4) 


335 


1125378 


651/2 


F 


— 


4>-> 


6^3/2 5^' 


= 3) 


351 


7215083 


65*1/2 


F 


= 


4>- 


6^3/2;^ 


= 4) 


351 


7217095 


65*1/2 


F 


= 


4>-> 


6^3/2 5^' 


= 5) 


351 


7219605 


6S1/2 


F 




4>-> 


rpx/2) 




652 


50476 


6S1/2 


F 


= 


4)- 


7P3/2) 




657 


93238 


65*1/2 


F 




4>-> 


8PI/2} 




770 


73660 


65i /2 


F 


= 


4>- 


8P3/2} 




773 


21409 


65*1/2 


F 


= 


3)- 


QPi/2-,F' 


= 3) 


335 


1205628 


65*1/2 


F 


= 


3)- 


QPi/2-,F' 


= 4) 


335 


1217305 


65*1/2 


F 


= 


3>- 


QP3/2-,F' 


= 2) 


351 


7305497 


65*1/2 


F 


= 


3)- 


QP 3 /2;F' 


= 3) 


351 


7307010 


65*1/2 


F 




3>- 


QPy2;F> 


= 4) 


351 


7309021 


65*1/2 


F 




3)- 


7P1/2) 




652 


51395 


65*1/2 


F 




3>- 


7P3/2) 




657 


94157 


65*1/2 


F 




3)- 


8P1/2) 




770 


74579 


65*1/2 


F 




3>- 


8^3/2) 




773 


22328 



These frequencies are known with an accuracy high enough so that we do not consider their 
contribution to the evaluation of the uncertainty in the final result. 

As far as the a n jF, b n jp and c n jF coefficients are concerned, we observe that in their definition 
the knowledge of the wave function at the origin is required. In turn, the wave function depends 
on the potential used to describe the interaction of the valence electron with the nucleus and with 
the closed electron shells |33| . In other words, this potential accounts for the departure from the 
hydrogen potential. Due to the difficulty in modelizing this interaction, we made a conservative 
estimation, attributing an uncertainty of 10% to the coefficients a n jF, b n jF and c u jf- 

Inserting all these quantities and the respective uncertainties in Q17|) we obtain: 



In (|27|) the uncertainty is mainly due to the uncertainty of the coefficients a n jF, b n jF and c n jF- 
We notice that the result (|27[) has been obtained by means of the integration in (|22|) over the 




(27) 
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solid angle 0, which has allowed us to calculate the scalar polarizability we are interested in; as 
a check, we obtain the same value using directly the Sandars formula |23| . If we estimate the 
BBR shift by using the polarizability value computed for a well defined orientation of the electric 
field with respect to the quantization axis, the non-zero tensorial part of the polarizability itself 
leads to small errors in the coefficient j3. More precisely, changes of about 3 x 10 -16 in 133 Cs, 
depending on whether the electric field is parallel or perpendicular to the z axis. 
^From (|17|) and 1)18|) we have also for the Stark shift coefficient: 

k = -(1.97 ± 0.09) x 10" 10 Hz/(V/m) 2 (28) 

and for the differential polarizability: 

aio = (1.14 ± 0.05) x 10~ 43 J/(V/m) 2 (29) 

Moreover, our computations give for the scalar polarizability «q of the 6 2 S , 1/ / 2 state of 133 Cs due 
only to the valence states: 

a v Q = (6.341 ± 0.016) x 10" 39 Cm 2 /V 

= (56.99 ± 0.14) x 10~ 24 cm 3 (30) 
= (384.9 ± 0.9)ajj 

In 1)30(1 the first result is the value expressed in SI units, the second one in CGS units and the 
third one in atomic units, being ao the Bohr radius. 

Finally, our evaluation of the corrective term e is (see Appendix) : 

e = 1.4 x 10~ 2 

in agreement with the value reported in 

IV. Discussion of the results 

We firstly examine the scalar polarizability value reported in ()30() : even if it does not lead di- 
rectly to the BBR shift of the hyperfine transition we are looking for, it represents anyway an 
important check of a part of our computations and provides an independent estimation of this 
important parameter. The total scalar polarizability ao must take the core polarizability a^ into 
consideration which accounts for approximately the 4 % of the total value in the case of 133 Cs 
[221 ■ O ur fi na l result for ao is then: 

a = a v + a c =(6.600 ± 0.016) x 10~ 39 Cm 2 /V 

=(59.32 ± 0.14) x 10 -24 cm 3 (31) 
=(400.7 ± 1.0)ajj 
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Table 6: Cesium ground state polarizability. 



ao( x 10~ 39 Cm 2 /V) Reference 

6.587 ± 0.031 Derevianko et al. [22] 

6.596 ± 0.013 Derevianko and Porsev |21| 

6.600 ± 0.016 this work 

6.611 ± 0.009 Amini et al. Q7j 



where is given by (|30|) and «q is reported in [30] . 

In Table 6 we report the data found in literature with an uncertainty lower than 1 %. 

An agreement at a level of ±1 x 10 -3 can be observed, well inside the quoted error bars, both 
between the theoretical computations and with respect to the experimental value of Amini and 
Gould PU- 

In Figure 2 the result of our computations concerning the BBR shift parameter (5 is compared 
with the experimental results reported in literature. 

Our result agrees with the direct measured values while the indirect measurements, even if 
in agreement among them, are higher by three standard deviations. This discrepancy has to 
be deeply examined not only for the correct estimation of the differential polarizability of the 
ground state hyperfine transition but also for its impact into the accuracy evaluation of the 
primary frequency standards. First of all, the following points have to be considered: 

(i) our theoretical result leads to a polarizability ao of the 2 5' 1 / 2 state in agreement with the 
most recent and precise evaluations as shown in Table 6, while the indirect measurements, 
considering that the Stark effect is proportional to the polarizability at the first order, lead 
to a considerably higher value for ao; 

(ii) our computations lead to the result of Itano et al. 31 if we use the old values of the transition 
frequencies [31] and of the electric dipole moments [201 ESI; the latter, in particular, are 
considerably higher than the most precise recent values; 

(iii) our (3 value agrees with the values obtained with the direct measurement method. 

Taking the above remarks into account, the following hypothesis may be considered in order 
to explain this discrepancy: 

(i) the assumption of equivalence between ac and dc Stark shifts is not fully correct, in spite 
also of the correction term e introduced by Itano et al. [3]; 
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This work 
a Levi et al. 



b Simon et al. 



c Bauch and Schroder 



d Mowat 



8 Haun and Zacharias 



i 1 1 1 1 1 1 1 

1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9x10" 



Figure 2: BBR shift (parameter (3) of the ground state 133 Cs hyperfine transition. With the 
triangle we indicated the indirect measurements, while with the square the direct ones. a is Ref. 
®, b is Ref. 0, c is Ref. 5,, d is Ref. 0, and c is Ref. 0. 



(ii) a systematic shift is present in the indirect measurements due to a physical effect not 
considered and present in the experiment. 

Since we have no reason at the moment to believe in the breakdown of the equivalence between 
ac and dc Stark shifts, as far as the second hypothesis is concerned, the following point should 
be examined. All the measurements performed with the indirect method are based on a Ramsey 
interaction scheme with a Stark field applied during the free flight of the atoms between the two 
interaction regions. The atoms experience then two sudden changes of the electric field which 
could perturb the free evolution of the hyperfine coherence or induce Majorana-type transitions. 

The last hypotheses must be considered more a clue than a probe to explain the discrepancy 
under question and obviously requires a more complete theoretical treatment. 

V. Conclusions 

We have reported in this paper a theoretical evaluation of the 2 £1/2 ground state polarizabil- 
ity of 133 Cs and the BBR shift of the ground state hyperfine transition between the magnetic 
field independent sublevels. The former is in really good agreement with recent theoretical and 
experimental values; the latter agrees with the direct measurements reported in the literature 
and is lower by three standard deviations than the values obtained from dc Stark measurements. 
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The discrepancy of 2 x 10 at T=300 K between the more recent estimations and the value 
accepted up to now impacts the accuracy evaluation of the primary frequency standards at the 
level of two or three times their believed standard uncertainties. Moreover, serious problems may 
be expected from the BBR shift to asses the accuracy at the level of 1 x 10~ 16 for the future 
microwave frequency standards operating at T ~ 300 K. 



Appendix 

To avoid the singularity at v = v ak in section II we have omitted the v 2 term in the denominator 
of (|7|). This simplification, usually adopted in literature, is justified by the fact that near room 
temperature the energy content of the BBR is peaked at frequencies much lower than the lowest 
transition frequency of the atom. An estimation of the error due to this approximation is provided 
in [Sj with the introduction of the corrective term e in ©. Here we re-evaluate this term following 
also the method described in [35] . 
Relations 0, © and © give: 

^c 3 V * J tLJo e'-l 
where ^ ak = |^ » 1 and x = -^p. Expanding in (|3*2*|) the factor: 

1 x 2 x 4 . , 

—-^1 + — + — + ... (33) 

1 ^ak ^ak 



and integrating we obtain: 

2 00 / , m\2n , 

„3 



Aw{k) _ 87r\d ak \ 2 u ak f k B T \ 2 ^ fk B T\ 2n (2n) 2n+2 \B 2n+2 
£oC 3 { h ) 2^\hu ak ) 4(n + l) 



(34) 



being B n the Bernoulli numbers. If we consider the main (n = 1) and the first order (n = 2) 
terms of (|34jl and sum over all the valence states, we obtain: 

We substitute now 1)35(1 in (|16j) and after some algebraic arrangements we obtain: 



4aT 4 fT\U^ \dp k \ 2 ^\d ak \ 2 



h 2 e c\T J vpk ^ Vak 

where 



, T 



2-1 



21 V h j r fel 2 _ K fc | 2 

Z^fe ^ L^k u ak 



(36) 



(37) 
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The main term of (|36|) coincides with Q17|) . as expected. As regards the corrective term e, following 
the same procedure to evaluate the dipole matrix elements in the basis of modified eigenfunctions 
as reported in section III, we obtain e = 1.4 x 10~ 2 , in agreement with the value reported in 
0. This agreement is not surprising because the recent values of the the dipole moments affect 
significantly the main parameter j3, but much less the corrective term e, as can be observed in 
(|37|) where they appear both in the numerator and in the denominator. 
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